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ABSTRACT 

The  computation  of  scattered  fields  from  fluid-loaded  bounded  objects  can  be  treated  in  a  consistent  manner 
using  the  extended  boundary  condition  (EBC)  or  T-matrix  method  for  electromagnetic,  elastic,  and  acoustic 
problems.  We  outline  the  general  principles  of  the  EBC  method  and  apply  it  to  problems  of  engineering  interest 
which  pertain  to  acoustical  scattering  from  submerged  three-dimensional  elongated  impenetrable  objects  or  elastic 
solids  and  shells.  Several  physical  examples  of  interest  are  offered  and  interpreted. 

INTRODUCTION 

Several  techniques  are  available  for  describing  waves  that  scatter  from  objects  of  known  constitution  and 
geometry.  Many  of  them  are  rather  specific  or  have  intractable  numerical  pitfalls.  It  is  therefore  desirable  to  obtain 
a  formulation  thar  allows  for  general  objects,  frequency  ranges,  and  boundary  conditions,  and  that  overcomes 
numerical  difficulties  commonly  encountered  in  several  broad  classes  of  numerical  methods.  In  this  paper  we  briefly 
describe  &  consistent,  unified  and  manageable  numerical  approach  useful  for  researchers  interested  in  solving  any 
of  a  wide  class  of  scattering  problems.  It  is  based  on  the  coupling  of  the  exterior  and  interior  solutions  of  the  surface 
boundary  representadon  of  the  Helmholtz  or  elastodynamic  equations  and  yields  the  Extended  Boundary  Condition 
(EBC)  method  of  Waterman1*1.  The  EBC  method  avoids  numerical  problems  often  encountered  by  other  techniques. 
We  present  numerous  physical  examples  which  are  chosen  not  only  for  their  intrinsic  interest  but  also  because  they 
represent  comparatively  difficult  problems  to  solve  by  other  means. 

EBC  FORMULATION 

The  EBC  technique  is  a  boundary  integral  method  that  couples  the  exterior  direct  scattering  field  solution 
'vith  the  intenor  field  solution.  The  method  was  initially  developed  for  electromagnetic  scattering  by  Waterman1  in 
1965,  and  subsequently  for  acoustical  scattering  by  him2  in  1969.  The  method  has  since  evolved  and  is  employed 
w  all  areas  of  classic  scattering  at  various  levels  of  sophistication5'1.  As  with  other  methods,  it  is  often  possible  to 
®ake  both  structural  and  numerical  variations  in  the  basic  EBC  approach.  We  have  developed  formulations5-5  which 
®°ur  experience  have  proven  the  most  favorable  for  numerical  implementation  and  we  employ  them  in  this  work. 
c  '*dth  the  simplest  formulation,  for  scattering  from  an  impenetrable  object  of  rotation. 

The  Huygen’s-Poincare’  integral  representation  of  the  scattered  field  U  exterior  to  the  bounded  object,  can 
oe  described  by  the  following  equation: 

-  Ut(r)-*-;§  fu+(r  )dG(r,r,y9n-G(r,r')3U+{rT/3n]dS.  (1) 
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The  quantity,  U.  is  the  scalar  wavefield  on  the  object  surface,  G  is  the  outgoing  Green's  function,  r'  is  taken 
on  the  surface  of  the  object,  and  n  is  a  unit  outward  normal  vector  to  the  surface  of  the  object.  S  is  taken  to  be  the 
surface  of  the  bounded  object.  In  the  subsequent  development  for  elastic  targets,  U.  will  be  the  displacement  vector, 
and  G  the  Green's  dyadic.  An  additional  expression  required  can  be  obtained  by  expressing  the  Huygen's-Poincare’ 
relation  at  a  point  interior  to  the  object,  which  leads  to  the  following  equation: 

0  - 11,(0+/.  [U,(0aG(rH.r')/3n-G(r".r')3Ut(0/an]dS  (2) 

where  r"  is  an  object  interior  point.  In  this  case,  the  field  can  be  described  as  being  "nulledM;  hence,  the  nomenclature 
"null  field"  method. 

Eqs.  (1)  and  (2)  yield  the  extended  boundary  condition  equations.  Since  in  their  present  form  they  are  not 
directly  useful,  we  now  reduce  them  to  a  form  amenable  to  numerical  computation.  Elastic  targets  submerged 
in  a  fluid  require  far  greater  mathematical  detail,  as  indicated  by  Pao‘  and  the  Varadans’.  Let  us  assume 
that  SU^frO/Sn  =  0,  so  that  we  obtain  the  expressions 

U(r)  =  U,(r)+(.  [Ut(iOaG(r,  r*)/9nldS  (3) 


and 


0  =  U,(r)  +/.  [U,(r')3G(r,  r'y9n]dS . 


(4) 


To  solve  these  expressions,  it  is  convenient  to  represent  U,(r),  UJr')  and  G(r.r’)  in  some  suitable  series 
expansion,  which  upon  truncation  leads  to  matrix  equations  that  can  be  solved  using  digital  computers.  The  Green  s 
function,  G,  is  a  normal  operator,  and  thus  can  be  represented  by  the  biorthogonal  series 


G(rj') « ikIRe  <p,(r<)<p,(r>) 


(5) 


where  r<and  r>  are  the  greater  and  lesser  of  the  two  points  r  and  r  ‘  relative  to  the  origin  of  the  object.  U,,  the  incident  j 
wavefield,  is  known  and  can  be  shown  to  be: 


U,(r)  =  la.Retp.fr). 

Similarly  the  scattered  field  U§(r)  =*  U(r)  -Ut(r)  can  be  expanded  as: 
U.(D  =  If.tp.(r) 
which  leads  to 

at=  -ik/U/ri)  3<p,(r0/3n  dS. 
and 


(6) 
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(7) 


f.“  ik(U.(ri)  9Retp,(ri)/3n  dS.  (T) 

The  following  expansion  has  been  proven  to  satisfy  closure  on  the  surface  for  the  rigid  problem: 

U.(r)  =  2.  b.Re<p.(r)  *  (8) 

where  bt  is  the  only  unknown  on  the  right  hand  side  of  Ea.  (7).  This  reduces  (7)  to  an  expression  in  which  tbej 
expansion  coefficients  b,  are  the  only  quantities  to  be  determined.  We  obtain 


*.  *  * b.  IVm  (r^  Rc<P.(0/9n  dS»  -Im  bm  Q^, 
f.=  2m  bm  /Re<p_  (r*)3  Retpi(rr)/3n  dS=  2m  b_  ReQ_, 


(9) 

(91) 
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where  is  in  element  of  some  known  matrix.  In  matrix  form.  Equations  (9)  and  (9*)  can  be  expressed 

a  =  -ikQb  (10> 

f  « ik  Re  Qb.  (ll> 

There  are  a  variety  of  numerical  strategies  to  attack  this  computational  problem.  The  most  straightforward 
(but  unfortunately  the  most  problematic)  method  would  be  direct  numerical  inversion,  to  eliminate  b  from  Eq.  ( 10) 
and  replace  it  in  Eq.  (1 1).  to  arrive  at  a  formally  convenient  expression  for  the  scattered  wavefield,  represented  by  f. 
This  leads  to  Eq.  (12),  relating  the  scattered  field  f  with  the  incident  field  a. 


f*-ReQQ-'a.  (12) 

Q  is  an  n  x  n  matrix,  and  a  and  f  are  l  x  n  column  vectors.  We  can  view  this  as  a  mapping  from  the  incident 
field  a  to  the  scattered  field  f  via  the  quantity  -RcOQ-1.  It  is  conventional  in  certain  areas  of  physics  to  describe  such 
a  mapping  as  a  transition;  hence  this  computational  mapping  is  often  referred  to  as  a  transition  or  "T  matrix'’ 
operation. 

As  noted,  a  straightforward  evaluation  of  Eq.  (12)  is  usually  not  the  most  efficient  way  to  solve  for  f.  An 
alternate  method  devised  by  us  to  obviate  this  difficulty  is  the  coupled  higher-order  T-matrix5-4  method  described 
below.  The  most  elementary  formulation  of  this  solution  is  to  block  Q  into  four  submatrices  (Q(ij),  where  ij  =  (1,  m) 
and  demand  that  the  Q(ii)  matrices  be  square.  Tliis  results  in  a  new  T  matrix  of  order  (1, 1)  which  is  convergent  and 
requires  an  inverse  no  larger  than  the  largest  pair  of  (ii). 

T(ll)  *=  -PG'1  (13) 


where 

P  -  Re  Q(U)  -  Re  Q(lm)  Q-'(mm)  Q(ml)  (14) 

and 


G  «  QOD-QClmJQYmm)  Q(ml).  (15) 

Eq.  (10)  is  related  to  an  integral  equation  of  the  first  kind,  and  can  be  transformed  using  a  method  developed 
by  us  into  an  eigenvalue  problem.  The  unknown  quantity  in  Eq.  (10)  is  the  surface  term  associated  with  the  column 
vector  b.  By  obtaining  b  directly,  we  then  have  the  solution  for  the  scattered  field  f  via  Eq.  ( 1 1 ).  To  accomplish  this 
we  create  a  Hermidan  matrix  H  «  Q  Qf,  where  Qf  is  the  adjoint  of  Q.  (It  is  easy  to  show  that  H  is  self-adjoint  or 
Hermidan).  The  eigenvalues  X,  and  the  eigenvectors  can  be  obtained  in  a  straightforward  manner,  described 
below,  where 

HP..X.P,.  (16) 

It  is  known  that  the  eigenvalues  \  of  Hermidan  operators  are  posidve,  real  and  monotonically  increasing, 
and  that  the  eigenvectors  p.  form  an  orthononnal  set  of  vectors.  Moreover,  it  can  be  shown  that  these  vectors  span 
the  space  associated  with  this  problem.  Thus,  we  can  expand  the  unknown  surface  quantity  b  in  terms  of  the  (5,'s. 
By  exploiting  the  orthonormality  of  the  eigenfunctions  of  H,  we  arrive  at  the  following  expression  for  f: 

f«-RtQP  lA(QP)f.  '  (17) 

The  more  complicated  case  of  elastic  targets  imbedded  in  a  fluid  may  also  be  derived,  in  which  many  of  the 
Previous  procedures  can  be  used.  Let  us  outline  the  solution  of  this  problem  for  an  elastic  solid.  We  still  require  the 
Huygen’s-Poincare’  integral  in  the  fluid  due  to  scattering  from  the  target.  However,  we  also  require  a  generalization 
°f  the  expression  for  an  object  in  an  elastic  environment.  An  expression  suitable  for  the  problem  we  wish  to  develop 
keen  obtained  by  Varatharajulu  and  Pao*  in  a  form  appropriate  to  derive  a  transition  matrix  for  elastic  targets, 
•dll  be  used  in  the  following  equations.  For  the  fluid  we  have; 
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U(r’)  =  Uj(r')  +  f[tJJ  -  t(U)  UJ  dS  (interior  point) 


(18) 


0  =  U(r")  +  /(tJJ  -  t(U)  UJ  dS  (exterior  point).  (19) 

Here  traction  t=n  •  (XI  VU+  p(VU+  UV)).  In  the  solid,  we  have  an  analogous  pair  of  equations 

U(r”)  =  U,(rO  +  /  [t_U  -  t(U)  UJ  dS  (interior  poin')  (20) 

0  =  U,(i0  +  J[t_U  -  t(U)  UJ  dS  (exterio.-  point)  (21) 

where  r"  is  in  the  object  interior  and  f  the  exterior.  t_  (tj  and  U_  (U  J  are  the  traction  and  displacement  vectors  on 
the  surface  as  approached  from  the  object's  interior  (exterior). 

Imposed  on  these  equations  we  have  the  following  boundary  conditions: 

U,.n  =  U_.n  (22) 

Vnsst-*n  (23) 

t  x  n  «  0.  (24) 

In  the  following,  we  consider  the  case  of  general  three-dimensional  solids,  using  a  spherical  coordinate 
representation.  Here  the  outward  wave  solutions  to  the  full  elastic  wave  equation  arc 

yl^  »  (K/ky*  l/k  V  (hj(kr)  (0,<t>))  (25) 

=  l/a0+l),/l  V  x  (rh,(kr)  (0,<J»)  (26) 

(27) 

* 

where  h,  arc  spherical  Hankel  functions  of  the  first  kind,  and  Y“'s  arc  spherical  harmonics  (normalized  to  unity). 
The  total  wavenumbers  of  longitudinal  and  transverse  waves  arc  k  and  k  Normalization  is  chosen  so  that  each 
wavefunction  is  normalized  to  a  unit  flux  for  any  closed  surface  that  encloses  the  origin.  For  a  nonviscous  fluid,  only 
k  is  needed,  so  that  (ic/k)  is  replaced  by  unity  in  Eq.  (26). 

As  in  the  previous  T-matrix  formulation  for  acoustic  targets,  the  general  technique  of  obtaining  a  transition- 
matrix  for  an  elastic  solid  is  to  eliminate  the  unknown  surface  wavefields  and  tractions.  This  elimination  is 
accomplished  by  judicious  use  of  the  boundary  conditions,  the  partial  wavefield  expansions  (Huygen's  principle  for 
waves  interior  to  the  scatterer).  The  original  derivation  of  the  following  expression  is  given  by  Bostrora*.  By  suitable 
elimination,  one  can  arrive  at  relation  (28),  between  the  incident  field  partial  wave  coefficients  { a )  and  scattered  field 
partial  wave  coefficients  (f): 

f  *  -QJ  R',P(QjR‘,P)',a.  (28) 

For  completeness  we  list  the  results  for  an  elastic  shell7: 

f  =  -Q*  P(Q^M''P)''a  *  (29) 

where  M  =  Q*+RT2+iTj,  and  Q,  P,  and  R  are  defined  by  the  following  surface  integrals: 


Q(;;)ij  *  k’/p,cuJ/  (X,V .  (Ou/Rejy,™  (n  •  (Ou/R e)y 
-  (n  •  (Ou/Re)]y,a,jn  •  t,  (Ou/Re)  y®}  •  dS 


(31) 


R(°)>j  =  kl/Pj®1/ t(Ou/Re)  Vi™]  •  [(Ou/Re)  y.^xn] 

-  { n  •  ((Ou/Re)  v„®K“  •  tj)((Ou/Re)  vyroxn] }  •  dS 


Py  *  k  JJl/P,a>1mn  *  ^  V  ®)n  •  V.0,l  •  dS  (32) 

V-ReQW  (33) 

Tj  corresponds  to  scattering  from  an  evacuated  inclusion  in  an  elastic  environment.  Without  going  into  detail 
regarding  the  solution  of  Eqs.  (28)  or  (29)  it  is  best  to  employ  the  unitary  method5. 

In  the  preceding,  we  have  outlined  methods  that  include  a  broad  spectrum  of  boundary  conditions  to  describe 
scattering  from  axially  symmetric  targets.  Implementing  these  expressions  is  not  without  pitfalls,  and  much  of  our 
past  efforts  has  focused  on  overcoming  them,  with  good  success.  These  methods  are  suitably  general  and  effective 
to  address  a  large  class  of  problems  without  encountering  complications  experienced  by  other  numerical  schemes. 

Below  we  apply  the  above  EBC/T-matrix  formulations  to  a  number  of  physically  illustrative  engineering 
problems,  that  underscore  some  of  the  more  dramatic  computational  features. 

APPLICATIONS  OF  THE  T-METHOD  TO  ENGINEERING  PROBLEMS 


We  first  treat  scattering  from  rigid  impenetrable  objects.  The  two  simplest  cases  are  for  spheroids  and  for 
cylinders  with  hemispherical  caps.  In  addition,  we  alsoconsider  the  case  of  scattering  from  elastic  solids  .where 
we  specifically  examine  resonance  phenomena. 


T-Mamx  Applications  to  Impenetrable  Problems 

There  are  two  classes  of  targets  for  impenetrable  problems,  Le.,  soft  and  hard  scattered.  They  do  not  support 
body  resonances;  therefore;  we  examine  acoustic  quantities  appropriate  for  nonresonant  targets  such  as  angular 
distributions  which  are  dependent  on  target  geometry  and  can  be  useful  to  determine  such  features  of  target  shape 
as  symmetry  or  elongation.  In  particular,  reflection,  diffraction,  and  generalized  Snell's  law  behavior  can  be 
observed  as  curved-surface  analogs  for  the  plane-layered  case. 

Bistatic  angular  distributions  correspond  to  measurement  of  a  scattered  field  at  any  point  in  space  for  some 
incident  wave  fixed  relative  to  some  source-object  orientation.  In  Fig.  1.  we  examine  a  rigid  spheroid  of  aspect 
(length-to-width)  ratio  of  16:1.  Ftg.  1  a-d  represent  scattering  from  the  object  along  the  axis  of  symmetry  (end-on) 
(a).  30  (b),  60  (c)  and  90  (d)  degrees  relative  to  the  symmetry  axis  (broadside).  The  values  of  the  incident  wavefield 
frequency  are  expressed  using  the  dimensionless  quantity  lcL/2,  where  L  is  the  object  length  and  k  the  total 
wavenumber  (k=2it/X.).  The  value  of  kL/2  in  Fig.  2  is  200,  which  implies  that  the  object  is  about  70  wavelengths 
long  and  thus  in  the  intermediate-frequency  region  where  neither  low  nor  high  frequency  approximations  apply.  In 
ah  figures  frequency  is  sufficiendy  high  that  wave  diffraction  effects  are  significant  in'  the  forward  scattering 
direction.  Perhaps  the  most  interesting  feature  of  the  four  plots  (Fig.  1  (a-d))  corresponds  to  a  reflection  at  the  (fairly 
flat)  side  of  the  object  for  scattering  angles  of  30  and  60  degrees.  This  reflection  can  occur  only  for  very  elongated 
objects  that  approach  flat  surfaces,  so  that  the  reflected  angle  is  almost  the  same  as  the  incident  ingle  (relative  to 
a  straight  line  through  the  axis  of  symmetry). 

EBC  Applications  to  Elastic  Solids 

We  now  examine  a  phenomenon  specific  to  elastic  objects  with  smooth  boundary  conditions  surrounded  by 
an  acoustic  fluid,  namely,  body  resonances.  The  body  resonances  examined  originate  from  the  curved-surface 
equivalents  of  seismic  interface  waves  of  pseudo-Rayleigh  or  Scholte  type,  propagating  circumferentially  to  form 
standing  waves  on  a  bounded  object  If  phase  velocities  are  slowly-varying  (as  a  function  of  frequency)  at  the  object 
surface,  resonances  occur  at  discrete  values  of  kL/2.  These  resonances  manifest  themselves  in  a  prescribed  manner 
(described  below).  For  elongated  elastic  solids,  three  distinct  resonance  types  occur.  The  first  kind  that  we  wish  to 
uiusoate  has  to  do  with  bending  modes  or  flexural  resonances.  For  unsupported  spheroids  a  plane  incident  wave  at 
5  degrees  relative  to  the  axis  of  symmetry  can  excite  these  modes  illustrated  in  Fig.  2a  and  2b  for  aspect  ratios  of 
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4  and  5  to  1.  It  can  be  shown*  that  the  lowest  mode  corresponds  to  2,  and  thereafter  3, 4,  etc.  The  interesting  aWI 
about  these  resonances  is  they  can  be  predicted  by  exact  bar  theories  and  coincide  nicely  with  result  herein 
particular  interest  is  the  effect  that  with  increasing  aspect  ratio  the  onset  of  resonances  occur  at  lower  lcL/2  val 
jl»  opposed  to  Rayleigh  resonances.  The  second  kind  (at  lower  frequencies)  are  due  to  leaky  Rayleigh  waves  andh°^N 
been  showTt  to  be  related  to  both  target  geometry  and  material  parameters  (notably  shear  modulus  and  dens' 
Resonances  can  in  this  case  best  be  observed  by  examining  the  back-scattered  echo  amplitude  and  phase  resnrml 
plotted  as  a  function  of  kL/2.  We  illustrate  this  for  an  aluminum  spheroid  of  aspect  ratio  10  to  1  at  end  on  incidenen'l 
Here  we  see  two  resonances  superimposed  on  the  semi  periodic  pattern  due  to  Franz  waves  associated  with  ti»M  t 
scattering.  If  we  subtract  rigid  scattering  (in  partial  wave  space)  from  the  elastic  response  then  we  are  left  with 
resonance  response  illustrated  in  Fig.  4  which  shows  that  for  aluminum  the  first  resonance  is  broad  while  the  secoodl 
diminishes,  in  contrast  to  WC  targets  which  yield  small  narrow  resonances  at  the  lower  value  and  a  larger  resonant 
for  the  next  value.  In  addition  to  the  above  wave  phenomena,  it  is  also  possible  to  excite  "whispering  gallery"  ^ 
resonances,  which  we  mention  next  Finally,  we  examine  broadside  resonances  for  a  4  to  1  steel  spheroid.  Here  we^! 
can  excite  three  phenomena.  At  the  lowest  value  we  can  see  a  spike  represenring  the  lowest  order  Rayleigh  resonance  ^ 
seen  end  on  corresponding  to  a  standing  wave  circumnavigating  the  largest  meridian  of  the  spheroid.  We  also  see 
weak  Franz  waves  similar  those  due  to  a  cylinder  and  than  we  see  the  lowest  order  Ra>  leigh  and  Whispering  Gallery  * 
resonances  corresponding  to  circumferential  waves  around  the  smallest  meridian. 
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Figure  1 .  Scattering  from  rigid  spheroid  with  aspect  ratio  of  16  to  lata  kU2  =  200  a)  end  on,b)30  degrees  relative 
to  axis  of  symmetry,  c)  60  degrees  relative  to  axis  of  symmetry  and  d)  broadside. 


2  4  «  8  10  14  ie  24  2  4  8  8  10  14  18  24 


Figure  2.  Flexural  resonance  predictedfrom  a  plane  wave  scattering  at  45  degrees  relative  to  the  axis  of  symmetry 
of  a  steel  spheroid  a)  of  aspect  ratio  4  and  b)  of  aspect  ratio  5. 


Figure  3.  Back  scattered  signal  from  aluminum  spheroid  of  aspect  ratio  of  10.  to  I  end  on  incidence. 


Figu/e  -i.  Back  scattered  signal  from  aluminum  spheroid  of  aspect  ratio  of  10  to  1  end  on  incidence  with  rigid 
background  subtracted. 


Figure  5.  Back  scattered  signal  from  steel  spheroid  of  aspect  ratio  of  4  to  l  broad  side  incidence. 
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